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Abstract 

Background: Medicago truncotulo Gaertn. (barrel medic) is cultivated as a pasture legume for its high protein content 
and ability to improve soils through nitrogen fixation. Toxic concentrations of the micronutrient Boron (B) in agricultural 
soils hamper the production of cereal and leguminous crops. In cereals, the genetic analysis of B tolerance has led to the 
development of molecular selection tools to introgress and maintain the B tolerance trait in breeding lines. There is a 
comparable need for selection tools in legumes that grow on these toxic soils, often in rotation with cereals. 

Results: Genetic variation for B tolerance in Medicago truncatula was utilised to generate two F 2 populations from 
crosses between tolerant and intolerant parents. Phenotyping under B stress revealed a close correlation between B 
tolerance and biomass production and a segregation ratio explained by a single dominant locus. M. truncatula 
homologues of the Arabidopsis major intrinsic protein (MIP) gene AtNIP5;1 and the efflux-type transporter gene AtBORl , 
both known for B transport, were identified and nearby molecular markers screened across F 2 lines to verify linkage 
with the B-tolerant phenotype. Most (95%) of the phenotypic variation could be explained by the SSR markers 
h2_6e22a and h2_21b19a, which flank a cluster of five predicted MIP genes on chromosome 4. Three CAPS markers 
{MtBtol-1,-2,-3) were developed to dissect the region further. Expression analysis of the five predicted MIPs indicated 
that only MtNIP3 was expressed when leaf tissue and roots were assessed. MtNIP3 showed low and equal expression in 
the roots of tolerant and intolerant lines but a 4-fold higher expression level in the leaves of B-tolerant cultivars. The 
expression profile correlates closely with the B concentration measured in the leaves and roots of tolerant and 
intolerant plants. Whereas no significant difference in B concentration exists between roots of tolerant and intolerant 
plants, the B concentration in the leaves of tolerant plants is less than half that of intolerant plants, which further 
supports MtNIP3 as the best candidate for the tolerance trait-defining gene in Medicago truncatula. 

Conclusion: The close linkage of the MtNIP3 locus to B toxicity tolerance provides a source of molecular selection tools 
to pasture breeding programs. The economical importance of the locus warrants further investigation of the individual 
members of the MIP gene cluster in other pasture and in grain legumes. 

Keywords: Medicago truncatula, Abiotic stress, Legume, Boron toxicity tolerance, Genetic analysis, MtNIP3, Gene cluster, 
Differential expression 



Background 

Medicago truncatula Gaertn. (barrel medic; Fabaceae), 
native to the Mediterranean region, is a self-regenerating 
annual pasture species grown on over 4.5 million hect- 
ares in the dryland cereal/livestock zones in southern 
Australia [1]. M. truncatula is adapted to a wide range 
of soils from sandy loams to clays and prefers neutral to 
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alkaline soils (pH (water) > 6.5). When grown in rotation 
with cereal crops their ability to fix atmospheric nitrogen 
contributes significantly to the nitrogen balance of fol- 
lowing crops. This lessens the need for inorganic nitro- 
gen fertilisers and hence increases economic returns to 
producers and benefits the environment [2,3]. Addition- 
ally as a non-host of many important cereal diseases, 
they are able to provide a valuable disease-break from 
pathogens such as Take-all {Gaeumannomyces graminis), 
crown rot (Fusarium graminearum) and cereal cyst 
nematode (Heterodera avenae) [4], Together with Lotus 
japonicus, M. truncatula is also a model organism for 
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legume biology with an international consortium provid- 
ing researchers with numerous genetic resources includ- 
ing two drafts of the sequenced genome and genetic 
maps (http://www.medicagohapmap.org). Its relatively 
small genome size (ca. 500 Mbp) and high levels of syn- 
teny to major crop legumes such as pea (Pisum sativum) 
and alfalfa (M. sativa) [5] make it an ideal system to 
conduct molecular genetics analyses on biotic and abi- 
otic stresses affecting legumes. 

A major problem facing M. truncatula and other le- 
gumes in arid world regions such as southern Australia, 
northern Africa, and the Mediterranean is the presence 
of phytotoxic levels of Boron (B) in the soil which can 
significantly impact on seed yield and quality [6]. Be- 
cause remediation of B-toxic soils is impractical [7], in- 
vestigators have looked to improve plant tolerance to B. 
In cereal crop species such as wheat and barley, genetic 
variation for B tolerance has led to the identification of 
quantitative trait loci (QTL) linked to B toxicity toler- 
ance [8], and more recently to the genes underlying 
these QTL [9-11]. This has facilitated the development 
of closely linked molecular markers that can be used in 
breeding programs to improve the tolerance of wheat 
and barley grown in areas where B toxicity is a problem. 
In M. truncatula, genetic variation for B toxicity toler- 
ance among different cultivars had earlier been reported 
[6] but so far has not been exploited to identify the gene 
(s) or genetic regions controlling this trait. 

Recent studies in barley and Arabidopsis offer clues as 
to what genes may control B toxicity tolerance in M. 
truncatula. These studies have implicated the involve- 
ment of genes encoding one of two types of proteins in- 
volved in B uptake and translocation: (1) efflux- type B 
transporters or (2) members of the major intrinsic pro- 
tein (MIP) family. The first efflux-type B transporter 
gene identified was the Arabidopsis AtBORl gene [12]. 
Although AtBORl conferred tolerance to plants under B 
deficient conditions, its homologue in barley, HvBotl 
(identical to HvBOR2), enabled plants to tolerate high 
levels of B [9,11]. The over-expression of an AtBORl 
paralog, AtBOR4, in transgenic Arabidopsis plants also 
increased their tolerance to high B levels [13]. AtNIP5;l 
was the first B transporter gene identified that encodes a 
member of the MIP family [14]. MIPs, also known as 
aquaporins, are channels for water and/or small non- 
charged molecules [15] and are comprised of four sub- 
families; the nodulin 26-like intrinsic proteins (MPs), 
the plasma membrane intrinsic proteins (PIPs), the tono- 
plast intrinsic proteins (TIPs), and the small basic intrin- 
sic proteins (SIPs). AtNIP5;l encodes a NIP subfamily 
member, and like AtBORl, was shown to be associated 
with B deficiency tolerance, with AtNIP5;l functioning 
as a boric acid channel required for the efficient uptake 
of B in roots [14]. In contrast, a recently discovered NIP 



gene in barley, HvNIP2;l, was found to be an important 
determinant of B toxicity tolerance [10], with the au- 
thors proposing that B toxicity tolerance is mediated by 
reduced expression of HvNIP2;l to limit B uptake as 
well as by increased expression of HvBotl to remove B 
from roots and sensitive tissues. A MIP subfamily mem- 
ber, AtTIP5;l was also recently described conferring tol- 
erance to B toxicity when over-expressed in transgenic 
Arabidopsis plants [16]. 

In this study, a candidate gene approach combined with 
population genetics and phenotypic analysis was used to 
determine the molecular basis of B tolerance in M. 
truncatula. We report the identification of a homologous 
AtNIP5;l gene, MtNIP3, in M. truncatula associated with 
B toxicity tolerance. A gene-based marker derived from 
MtNIP3, named MtBtol-1, and further markers co- 
segregated with B toxicity tolerance in two F 2 populations. 
In addition, gene expression of MtNIP3 was up-regulated 
in the leaves of respective B-tolerant parental lines. This 
information warrants further characterisation of MtNIP3 
as a candidate gene for B tolerance and provides gene- 
based markers as molecular tools to introgress and main- 
tain the tolerance trait in Medicago breeding programs. 

Results 

Segregation analysis 

The inheritance results of the phenotypic data are 
presented in Table 1. The segregation ratio for the B tox- 
icity tolerance trait has a significant fit (p = 0.05) to the 
expected 3:1 (tolerant vs. intolerant) ratio for a single 
dominant locus, which is true for both F 2 populations. 
F 2 plants of the population tap x Paraggio (intolerant x 
tolerant) had leaf symptom scores ranging from 0 to 4, 
whereas in tap x Caliph the scores ranged from 0 to 5. 
Relative to these symptom ranges in txP and £xC, plants 
with a score of 0 and 1 in txP and a score of 0 to 2 in 
txC were considered as B tolerant. 

B accumulation in root and leaf tissue 

Initially, subsets of 39 and 36 F 2 lines from the txP and 
txC populations, respectively, were selected to screen for 
potential molecular markers. The subsets consisted of 
lines that were phenotyped after three weeks of high B 
treatment and had clear tolerant (0 or 1) or intolerant 
(4 or 5) scores. ICP analysis was used to (a) quantify the 
amount of B accumulating in the leaf, and (b) verify the 
correlation between B levels and phenotypic scores. 
Additional ICP analysis was carried out on leaves and 
roots of the three parental lines to quantify the B distri- 
bution in tolerant and intolerant plants. The results are 
summarised in Figure 1. In the txP subset, the mean B 
concentration in leaves sampled from B-tolerant and - 
intolerant F 2 lines was 559 mg/kg (range: 420-710) and 
942 mg/kg (range: 760-1360), respectively. In the txC 
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Table 1 Phenotypic segregation of the B toxicity tolerance trait in the F 2 populations tap x Paraggio and tap x Caliph 



Population 


No. of plants phenotyped 


No. of B-tolerant F 2 plants 


No. of B-intolerant F 2 plants 


Segregation ratio x 2 


top x Paraggio 


316 


247 


69 


3 : 1 0.882 


top x Caliph 


313 


228 


85 


3 : 1 0.406 



subset the mean B concentration was 605 mg/kg (range: 
195-860) and 1021 mg/kg (range: 870-1210). The leaves 
of parental lines had mean B concentrations of 536 
(Paraggio), 592 (Caliph), and 1319 (tap) mg/kg. In con- 
trast to the mean B concentrations in the leaves, with 
tolerant lines accumulating less than half the amount of 
B compared to the intolerant lines, mean B concentra- 
tions in the roots did not differ significantly, 189 
(Paraggio), 201 (Caliph) and 200 (tap) mg/kg. 

Biomass production and Boron tolerance 

Breeding for B tolerance aims at maintaining yield under 
high B concentrations. As a proxy for the agronomic 
value of the B tolerance trait in pasture legumes, plant 
biomass or dry matter production (DM) was measured. 
The average DM production of the tolerant F 2 plants 
and their tolerant parent, Caliph, was not significantly 
different, 4.0 g versus 3.8 g, respectively. In contrast, the 
average DM of the intolerant F 2 lines was significantly 
higher than that of the intolerant parent tap, 1.9 g versus 
0.9 g, respectively (Figure 2A). The DM and symptom 
scores were highly correlated with a coefficient of -0.995 
(Figure 2B) indicating that the observed differences in 
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Figure 1 Distribution of B concentrations in tolerant and 
intolerant M.t lines after 3 weeks of excess B. The amount of B 
was quantified in leaves of tolerant and intolerant F 2 lines of both 
populations and the two tolerant parental lines Paraggio and Caliph 
and the intolerant parent top. ICP data were generated on 35-40 F 2 
lines of the tolerant and intolerant subsets of both populations, 4-8 
repeats of leaf and 4 repeats of root samples for each parent. 



biomass production are due to the B tolerance trait be- 
ing present or absent in the individual F 2 plants. With 
each increase in visual symptom score DM was reduced 
by approximately one gram. 

Development of molecular markers for B toxicity 
tolerance 

The AtBORl and AtNIP5;l cDNA sequences (GenBank: 
AB073713 and NM_117106) and their respective predicted 
protein sequences (GenBank: BAC20173 and NP_192776) 
were used to query the M. truncatula pseudomolecule 
database (version Mt3.5) available at the M. truncatula 
Hapmap project website (http://www.medicagohapmap. 
org). The best AtBORl homologue, identified on chromo- 
some 7 (BAC AC157892; GenelD: Medtr7gl 100001) en- 
codes a predicted protein, which shares 78% amino acid 
(aa) sequence identity with AtBORl. We developed a gene 
based CAPS marker for this gene but found no linkage be- 
tween this Boron transporter gene and B tolerance in ei- 
ther of the £xP or txC F 2 subsets. 

Four Medicago database matches to AtNIP5;l were iden- 
tified on chromosomes 1 (BACs AC181971 and AC225543) 
and 4 (BACs AC168305 and AC153003). Both chromo- 
some 1 candidates (GenelD: Medtrlg0978401 and 
Medtrlg0980801) encode exactly the same predicted pro- 
tein, which shares 69% sequence identity with AtNIP5;l, 
whereas the chromosome 4 candidates encode different 
predicted proteins with 60% (GenelD: Medtr4g0067301) 
and 54% (GenelD: Medtr4g0066301) identity to AtNIP5;l, 
respectively. The Medtr4g0067301 protein product corre- 
sponds to the aquaporin MtNIP3 (GenBank: AY539749) 
identified [17] and its corresponding gene is referred to as 
MtNIP3. The two BACs on each chromosome were adja- 
cent to one another and the physical distance between can- 
didate genes was small - approximately 84 kb and 71 kb 
for chromosome 1 and 4 candidates, respectively. Two 
flanking SSR markers (002G02, h2_13clla) located within 
500 kb of the chromosome 1 candidate genes showed poly- 
morphisms between tap and Paraggio but proved not to be 
linked to B tolerance in the subset of F 2 lines tested. Thus, 
like the AtBORl homologue on chromosome 7, the 
chromosome 1 candidate genes were ruled out from further 
analysis. In contrast, SSR markers flanking the NIP candi- 
date genes on chromosome 4 showed strong linkage with 
the B tolerance phenotype. The SSR markers MTIC348, 
h2_13f22h, h2_6e22a, h2_21bl9a and 001D02, which flank 
a 6205 kb interval were polymorphic in either one or both 
F 2 populations and explained 94% to 96% of the phenotypic 
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Figure 2 Mean plant dry matter (DM) comparison. (A) DM of 

tolerant and intolerant F 2 lines is compared to tolerant and intolerant 

parental lines (Caliph, tap) of the corresponding population. Error bars 

represent standard deviations of the mean. (B) High negative 

correlation (-0.995) between DM and severity of the visual B stress 

symptom scores, n = 21, 42, 32, 30, 13 for plants with visual symptom 

score 1 , 2, 3, 4, 5, respectively. 
\ ) 

variation (Table 2, Figure 3). Amongst these SSRs, the 
markers which were more distant from the 71 kb NIP gene 
cluster were less closely linked to the tolerance phenotype 
within the segregating populations (Table 2, Figure 3). To- 
gether with the segregation ratio characteristic of a single 
dominant locus, the strong linkage of these SSR makers to 
the tolerant phenotypes indicated that the trait-defining 
gene(s) reside at this 71 kb NIP gene locus. 

To dissect this locus, we derived CAPS markers from or 
closer to the NIP genes themselves. The M. truncatula 
Hapmap project gene annotations revealed that the candi- 
date AtNIP5;l homologue on BAC AC153003 has a poten- 
tial splice variant, with both gene products (GenelD: 
Medtr4g006730.1 and Medtr4g0067302) predicted to 
encode MPs. Likewise, the second chromosome 4 candidate 
on BAC AC168305 (GenelD: Medtr4g006630A) is also 
predicted to encode a NIP, with the 71 kb genomic region 
in between the two genes (Medtr4g006730A and 
Medtr4g006630.1) containing three more potential NIP 
genes (GenelDs: Medtr4g006650 A/660 A/670 A). Therefore, 
according to the gene annotations, the locus as a whole con- 
tains up to five predicted NIP genes encoding six gene prod- 
ucts potentially involved in controlling B toxicity tolerance. 



Three CAPS markers, named MtBtol-1, MtBtol-2 and 
MtBtol-3 (Table 2, Figure 3, Figure 4A-C) were developed 
that were polymorphic in both txP (200 F 2 lines tested) and 
txC (36 F 2 lines tested) populations (Figure 4D, E) within 
this region and explained 96% and 94% of the phenotypic 
variation, respectively. MtBtol-1 is specific to the first candi- 
date AtNIP5;l homologue, Medtr4g006730 A/730.2, MtBtol- 
2 is specific to the Medtr4g006650A gene which is situated 
13 kb away from the second candidate AtNIP5;l 
homologue, Medtr4g006630A, and MtBtol-3 is specific to 
the Medtr4g006660A gene. It was not possible to design a 
gene specific marker for Medtr4g006630A due to its close 
similarity with the Medtr4g006730 A/730.2 gene. 

NIP gene expression 

We used Q-PCR to investigate the expression profiles of 
the NIP genes in roots and shoots of the three parental 
lines over a 48 h period following excess B treatment. Due 
to the high level of coding sequence similarity between 
the Medtr4g006630A/730A/730.2 gene products, it was 
only possible to design gene specific primers for 
Medtr4g006630A and the splice variant Medtr4g006730.2. 
Design of gene specific primers was also hindered by the 
absence of information related to the untranslated regions 
of these gene sequences. Hence, a third primer pair, spe- 
cific to all three genes, was designed to generate a cumula- 
tive expression profile of B transporter candidates. Each of 
the primer pairs amplified the expected sized fragment 
when genomic DNA was used as template, however, only 
the primer pair specific for all three gene products ampli- 
fied clear fragments of expected size from the cDNAs of 
the time course experiment. Two Medtr4g006630A spe- 
cific primer pairs with different annealing positions in the 
gene amplified no product at all whereas the 
Medtr4g006730.2 specific primers yielded a small amount 
of expected product corresponding to <10 copies/ uL in 
the leaf and root tissue examined, which is at the very 
limit of detection for Q-PCR. Therefore, by exclusion, the 
cumulative expression profile that was generated for all 
three genes was in fact presumed to represent one, 
Medtr4g006730A. 

The expression profile for Medtr40006730A showed 
that the gene is constitutively expressed in the leaves, 
with approximately 4 times more transcripts present in 
both B-tolerant lines compared to the intolerant line 
(Figure 5A). In roots, the expression was equally low 
and non-differential in all three genotypes irrespective of 
their B tolerance phenotype. The expression level in the 
roots was comparable with the leaf expression of the in- 
tolerant line tap (Figure 5B). The stronger gene expres- 
sion of Medtr40006730A in the leaves of the tolerant 
lines correlated with the B toxicity tolerance phenotype. 

Q-PCR assays were also designed to investigate the 
possible involvement of the other three predicted MIP 
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Table 2 Molecular markers and the phenotypic variation they explain in the two F 2 populations studied 



Marker Marker, Mt Primer sequence (forward and reverse, 5' - 3') Restriction Expected CAPS Phenotypic 

chromosome enzyme fragment variation 

size [bp] explained [%] 











tap 


Par/Cal 


fx P* 


fxC** 


MTIC348 


SSR, chr. 4 


TGGAGGAGGGGTAGGATAGG ATGATGATGAGGCGGAGAAG 








95 


94 


h2_13f22h 


SSR, chr. 4 


GGOTCCTGATGCTGGTOG ACAAGCAGGTOGACACACA 








95 


94 


h2_6e22a 


SSR, chr. 4 


CCATGGTGTCAmGATCCTCCATGGTGTCAmGATCCT 








96 


n.p. 


MtBtol-1 


CAPS (gene), chr. 4 


CTCmACCTCCCCCTCCTGTOAACTACGATCAGGGTCGTOCTA 


Apol 


387 


170/217 


96 


94 


MtBtol-2 


CAPS (gene), chr. 4 


toaactacgatcagggtcgtoctatgctccagcacatccaa™ 


Drol 


102/131 


102/256 


96 


94 










256/280 


411 






MtBtol-3 


CAPS (gene), chr. 4 


TCCCACTGAAAATCGATOAA GGCCAAAGCTCATOCA^ 


Mspl 


125/450 


575 


96 


94 


h2_21b19a 


SSR, chr. 4 


GGAACACCGGACACAATACC CAC1TCTCAGGCACAGCAAA 








96 


n.p. 


001 D02 


SSR, chr. 4 


TCCCAACGCI 1 1 1 ICAmC GAAOTGAAGAAGAACGCCG 








94 


89 


MtChr7 


CAPS, chr. 7 


TOGAGCAmATACCAGCAA TOTATGCATCGGAGACTGC 


Nlolll 


306 


84/222 


n.a. 


n.a. 


002G04 


SSR, chr. 1 


agggtcgcctcaacta™ TCAACACCA^CTCAATG 








n.a. 


n.a. 


002G02 


SSR, chr. 1 


TAACGTOCTCCCTCCTCCGTCTCCAACGAAGTOAAGGG 








n.a. 


n.a. 


h2_13c11a 


SSR, chr. 1 


TCACAACOTOTGCACCAGTCATAAOTGACCTGCACCA 








n.a. 


n.a. 



Most closely linked markers are in bold font. * n = 200; ** n = 36; n.a., not applicable; n.p., not polymorphic. 



genes at this locus, Medtr4g006650. 1/660, 1/670,1. Again, 
the specificity of the primers was confirmed with the 
correct sized PCR fragment being amplified from DNA 
template. However, none of the primer pairs amplified 
any product on the cDNAs of the time course experi- 
ment from leaf and root tissue, suggesting no detectable 
expression before or after B treatment. 

Analysis of MtNIP3 (Medtr4g006730A) gene copy number 

We used Southern hybridisation to determine if the 
presence of additional MtNIP3 gene copies in the two 
B-tolerant lines could account for the increased mRNA 
transcript levels observed for this gene in the leaves. As it 
was not possible to design a specific hybridisation probe 
for Medtr4g006730.1 (MtNIP3) due to the high level of 
homology it shares with Medtr4g0066301 {MtNIPS), the 
designed probe would detect both genes. The resulting hy- 
bridisation revealed the presence of two gene copies in 
tap (as expected) and Paraggio, and one copy in Caliph 
(Figure 6). Thus, the observed increase in MtNIP3 tran- 
script levels in the leaves of the two B-tolerant lines was 
not associated with a higher gene copy number. 

MtNIP3 coding sequence and promoter region 

The MtNIP3 coding region was sequenced from Caliph, 
Paraggio, and tap to identify potential differences in the 
predicted amino acid composition of the corresponding 
MIP protein. A 918 bp cDNA fragment was amplified 
from all three genotypes. None of the seven single nu- 
cleotide polymorphisms (SNPs) identified were specific 
to Caliph/ Paraggio or tap (Additional file 1). The seven 
SNPs translated into three changes among the 305 



M.t. chr. 4 
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Figure 3 Chromosome region linked to B tolerance. Scheme of 
the genome region on linkage group 4 in M. truncotulo (version 
Mt3.5) showing the physical position of the B tolerance candidate 
genes and markers. 
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Figure 4 CAPS markers linked to B toxicity tolerance. (A) MtBtol-1; (B) MtBtol-2; (C) MtBtol-3; Lane L - DNA ladder, Lane 1 - undigested 
amplicon from tap, Lanes 2 to 4 - digested amplicons from tap, Caliph, and Paraggio, respectively. (D) and (E), Application of MtBtol-1 on 
genomic DNA from 22 randomly selected F 2 plants derived from the tap x Paraggio and tap x Caliph crosses, respectively, distinguishing 
homozygous and heterozygous F 2 plants. 



amino acid residues of the NIP protein in the three ge- 
notypes. Again, the changes were not specific to the 
B- tolerant genotypes (Additional file 2). 

A comparative Medtr4g0067301 promoter analysis be- 
tween Caliph, Paraggio, and tap was also performed to 
search for differences which could explain the variation in 
transcript levels observed between the B-tolerant and - 
intolerant lines. Sequencing of the 1.1 kb region directly up- 
stream of the transcription start site revealed a 99% basepair 
identity between the three genotypes. The differences were 
not specific to the intolerant or the tolerant genotypes and 
further work would be required to address the impact of 
these single and dinucleotide differences on gene expression 
or search for potential differences further upstream. 

Discussion 

In this study, a candidate gene approach was used to identify 
the underlying genetic loci for tolerance to B toxicity in M. 
truncatula. Similar to the recently reported identification of 
the genetic basis for herbicide tolerance in M. truncatula 
[18], we utilised the molecular and physiological information 
on B homeostasis first reported in Arabidopsis [12,14]. Al- 
though the corresponding proteins of the Arabidopsis B 
transporter /channel genes AtBORl and AtNIPS;l are asso- 
ciated with the efficient transport of B across the plasma 
membrane under B deficiency [12,14], similar proteins in 
barley were shown to aid the plant in conditions of B excess 
[9-11]. Therefore, it was hypothesised that homologous 
genes in M. truncatula of AtNIPS;! and AtBORl could also 
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Figure 5 Quantitative RT-PCR of MtNIP3. (A) MtNIP3 expression in 
Paraggio (P), Caliph (C), and tap (t) is regulated in the leaves (B) but 
not in roots following excess B treatment. Error bars represent 
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Figure 6 Southern hybridisation with an MtNIP3 
{Medtr4g006730A)-spec\f\c probe in Caliph (C), Paraggio (P) and 
tap (f). 



play a role in supporting plants to tolerate high levels of B in 
the soil 

Genetic variation for B toxicity tolerance was identi- 
fied in the M. truncatula lines Paraggio and Caliph (tol- 
erant, [6]) and tap (intolerant). Paraggio and Caliph 
were each crossed with tap to generate two F 2 popula- 
tions segregating for B tolerance. The phenotypic segre- 
gation of the B tolerance trait was measured in over 300 
individuals for each of the F 2 populations and agreed 
with a single dominant locus that might host a single or 
several closely linked tolerance genes. A very close cor- 
relation between the visual symptom scores on individ- 
ual plants and their biomass production was also evident 
suggesting B uptake as the direct and major cause of re- 
duced plant growth and validating the symptom scores 
as a valid tool to estimate plant performance. All toler- 
ant plants analysed showed a reduced (on average ca. 
50%) level of B in their leaf tissue compared to intolerant 
lines implying an exclusion mechanism being the under- 
lying B tolerance mechanism in this pasture legume. An 
exclusion mechanism was also proposed in barley where 
B-tolerant cultivars showed decreased B concentrations 
in their leaves compared to intolerant cultivars [11,19]. 

With the aid of the 2nd draft (version Mt3.5) of the se- 
quenced M. truncatula genome, candidate AtNIP5;l and 
AtBORl homologues were identified on chromosomes 1, 
4, and 7. A subset of highly B-tolerant and -intolerant 
lines from the two segregating F 2 populations was used 
initially for genotyping to establish linkage between the 
B-tolerant phenotype and SSR markers located close to 



the homologous Bor and NIP genes on all three chromo- 
somes. Only SSR markers located near two candidate 
AtNIP5;l homologues, separated by only 71 kb, on 
chromosome 4 showed linkages of approximately 95% to 
the tolerance trait in both F 2 populations. The strong 
linkage was then upheld across 200 individual F 2 lines 
from the tap x Paraggio cross which further supported 
that the corresponding genetic region was associated 
with B toxicity tolerance. We referred to one of these 
candidate genes as MtNIP3 and the predicted protein as 
MtNIP3, which was one of six aquaporins first identified 
in M. truncatula [17]. The other predicted gene was re- 
ferred to as MtNIPS. 

Although both corresponding proteins encoded by 
MtNIP3 and MtNIPS confirmed good amino acid se- 
quence similarity to AtNIP5;l, they both showed a 
greater level of amino acid sequence identity to 
AtNIP6;l - another MIP family member and B channel 
identified in Arabidopsis [20]. MtNIP3 has 75% se- 
quence identity to AtNIP6;l whereas MtNIP5 shares 
only 57% identity. Within the 71 kb interval, between 
MtNIP3 and MtNIPS, there are up to three additional 
predicted MIP genes and a splice variant of MtNIP3. 
Fine mapping this locus with three CAPS markers de- 
rived from the actual MIP genes themselves or se- 
quences nearby did not improve the linkage obtained 
with the two most closely linked SSR markers. However, 
the slight decrease in linkage for flanking SSR markers 
either side of this locus gave a strong indication that one 
or a combination of the MIP genes in the 71 kb interval 
were involved in regulating the transport of B in the 
plant. Therefore, as an output of the genetic analysis we 
identified eight molecular markers that were closely 
linked to the B tolerance trait, the best of which were 
CAPS markers MtBtol-1, MtBtol-2, and MtBtol-3 and 
SSR markers h2_6e22a and h2_21bl9a, which are all 
95% linked. The observation that these markers explain 
95% and not 100% of the phenotypic variation is 
expected considering that B-intolerant F 2 plants still 
show a greater level of B tolerance, when quantified as 
DM production, than the intolerant parent (Figure 2A). 
Plant characteristics independent of B exclusion such as 
morphological differences e.g. a stronger vigour of Ca- 
liph and Paraggio over the intolerant tap mutant, are 
likely to account for the remaining 5% of phenotypic 
variation in these two populations. Aside from DM pro- 
duction, the ICP analysis clearly showed differences in 
the range of B levels in the leaves of the F 2 lines from 
the txC compared to the txP population. These differ- 
ences can again be explained by genotype-specific traits 
(e.g. maturity, morphology) that are inherited independ- 
ently of the B exclusion locus. 

As the chromosome 4 locus associated with MtBtol-1 
contains up to five predicted MIP genes, we tried to 
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narrow down the association with B toxicity tolerance to 
individual genes by monitoring their expression in Ca- 
liph, Paraggio, and tap prior to and after treatment with 
excess B. Although a specific Q-PCR assay could not be 
designed for MtNIP3, by excluding the low or non- 
expression of its predicted splice variant and the similar 
gene MtNIPS, respectively, we deduced that MtNIP3 was 
the only MIP gene expressed at this locus as no expres- 
sion signals were detected for the remaining three MIP 
candidates when using alternative sets of primers. This 
meant that they were either not expressed, expressed at 
extremely low levels, or that the gene predictions were 
not correct and these genes were simply not there. Ul- 
timately, the presence and function of in silico predicted 
genes has to be confirmed through the analysis of their 
gene product [21]. Most interesting was the observation 
that the constitutive expression of MtNIP3 in the leaves 
was almost four-fold higher in both B-tolerant lines com- 
pared to the intolerant line. By contrast, there was no dif- 
ference in expression levels in the roots of tolerant and 
intolerant lines. A similar pattern was observed with the 
distribution of B in the leaves and roots of the tolerant 
and intolerant parental lines. There was no significant dif- 
ference of B in the roots but twice the B concentration in 
the leaves of the intolerant parent compared to that in 
both tolerant lines. Several mechanisms of alleviating B 
toxicity are being discussed, including efflux from the 
roots, redistribution of B in the leaves and leaching of B 
by rain or removing B from the plant through guttation 
[22]. Based on the observed B distribution in the roots 
and leaves we can exclude that efflux from the roots is re- 
sponsible for the different levels of B tolerance in these ge- 
notypes. Nevertheless, other B tolerant medic genotypes 
[6] might well employ such a mechanism. Possible mecha- 
nisms could be either the redistribution of B in the leaves 
from sensitive symplastic to apoplastic spaces in the leaves 
and subsequent leaching through rain [22] and/or the re- 
moval of B via guttation [9,22,23]. After B treatment, 
MtNIP3 expression was regulated in the leaves only, with 
the level of differential expression being 4 times higher in 
B-tolerant compared to intolerant lines. The fact that the 
two B-tolerant lines also have substantially less B in their 
leaves suggests that MtNIP3 itself or homologs at the 
MtNIP3 locus may function in removing B from the leaf. 

Southern hybridisation with an MtNIP3l MtNIPS specific 
probe ruled out the possibility that additional copies of 
MtNIP3 were responsible for the increased expression of 
this gene in the leaves of B-tolerant lines. In fact, Caliph 
only had one gene copy as opposed to the two detected in 
Paraggio and tap. In barley, it was shown that a B-tolerant 
genotype had approximately four times as many copies of 
the B transporter gene Botl compared to a number of in- 
tolerant genotypes [9], however, in Medicago this does not 
appear to be the case for MtNIP3. 



With MtNIPS not showing measurable transcript levels 
in roots and youngest leaves, our findings suggest that 
MtNIP3 is the most promising candidate to be involved in 
the mechanism controlling B toxicity tolerance in M. 
truncatula as its expression is strongly linked with B toler- 
ance and distribution in the plant. Given that B is essential 
for cell wall elongation, where it is required for the cross- 
linking of the pectic polysaccharide rhamnogalacturonan II 
[24], the constitutive expression of MtNIP3 in the roots 
and leaves is not surprising as MtNIP3 may always be 
needed to give B access to sink tissues for normal plant 
growth. The lower constitutive and then non-differential 
MtNIP3 expression observed in the roots following B treat- 
ment is consistent with what was reported for AtNIP6;l } 
the Arabidopsis gene most similar to MtNIP3 [20]. It 
would be interesting to know if this closest Medicago 
match to AtNIP6;l does share the water impermeability 
of the AtNIP6;l protein [20,25]. Considering that M. 
truncatula grows in low rainfall environments of northern 
Africa, Mediterranean Europe and southern Australia, a 
water conserving mechanism to remove excess B would 
be advantageous to the plant. It is likely that additional 
MtNIP3 alleles, homologous or orthologous genes to the 
efflux transporter genes AtBORl [12] or HvBotl [9,11] 
control B toxicity tolerance in other M. truncatula lines 
than the ones investigated here. The information gained in 
M. truncatula could also be used to facilitate the identifi- 
cation of genes in other pasture legumes (e.g. other 
Medicago species [6]) and also grain legumes (lentils [26], 
peas [6]) as recently demonstrated for alfalfa and the re- 
sistance to the fungal disease anthracnose [27]. 

The generated molecular markers are currently being used 
in SARDIs pasture breeding program to confirm that B tol- 
erance is conserved in advanced lines that also contain toler- 
ance to the sulfonylurea family of herbicides [28], which is 
being selected for by using the gene-specific marker CAng 
[18]. These breeding aims support a crop/pasture rotation 
system where self-regenerating annual medics are grown in 
rotation with wheat and barley to fix atmospheric nitrogen 
to lower input costs, provide a cereal pathogen disease- 
break and serve as a high protein feed for livestock. At an 
estimated rate of approximately 25 kg nitrogen per tonne 
DM and a production of 4 tonnes DM per hectare, a medic 
crop can fix approximately 100 kg nitrogen per hectare [29], 
reducing the need for expensive inorganic fertilisers and at 
the same time avoiding C0 2 emissions that would arise 
from the methane-driven fertiliser production. 

Conclusions 

In the present study, a candidate gene approach was ap- 
plied to establish a genetic linkage to the B toxicity toler- 
ance trait observed in different cultivars of Medicago 
truncatula. Whereas no genetic linkage between other 
Boron exclusion and efflux transporter genes was 
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detected, strong linkages were evident to genetic markers 
flanking a cluster comprising five predicted NIP genes 
explaining most (95%) of the phenotypic variation ob- 
served in the segregating Medicago populations. The de- 
rived markers are valuable tools for selection of this trait 
in a pasture breeding program and the close correlation of 
biomass and B toxicity qualifies B toxicity tolerance as a 
priority trait. The close correlation of MtNIP3 expression 
and B concentration in roots and leaves of tolerant and in- 
tolerant lines suggests that this member of the gene clus- 
ter is the most promising candidate for a B tolerance gene 
in this species. It further suggests that the molecular basis 
of B detoxification in Medicago differs from that observed 
in the monocot barley where differential expression of 
HvNIP2;l in the roots [10] and a greater Boron-efflux 
transporter gene activity in the leaves [9] has previously 
been identified as the molecular basis of B toxicity allevi- 
ation. Further physiological studies on guttation and B 
leaching through rain, as recently described for barley [22] 
are required to quantify the possible role of these two 
mechanisms in mediating B toxicity tolerance. 

Methods 

Genetic material 

Two F 2 populations of M. truncatula were developed by 
crossing two B toxicity tolerant cultivars, Caliph and 
Paraggio (Caliph tested as Z-602 [6]), with the B toxicity 
intolerant line, tap. Paraggio and Caliph are cultivars de- 
veloped by the South Australian Research and Develop- 
ment Institute (SARDI) and were registered as cultivars in 
1982 and 1993, respectively, whereas tap {mtapetala) is a 
male sterile mutant derived from EMS treatment of A17 
(Jemalong A17), the genotype used for M. truncatula gen- 
ome sequencing [30]. For determination of the B toler- 
ance trait inheritance, 316 and 313 F 2 plants were 
assessed from the tap x Paraggio (£xP) and tap x Caliph 
(txC) crosses, respectively, whereas 200 and 36 F 2 plants 
from the £xP and txC crosses, respectively, were used 
for molecular marker analysis. 

Plant growth conditions 

Seeds from the three parental lines and the F 2 progeny 
were surface sterilised, scarified with a small nick to the 
seedcoat with a scalpel blade, and allowed to germinate 
on moist filter paper in an inverted position. When the 
radicals had reached approximately 10 mm in length 
seeds were transferred onto a wire mesh and suspended 
over a 32 x 45 cm tub containing 10 L of a modified ver- 
sion of Hoaglands nutrient solution (2 mM KN0 3 , 2 
mM Ca(N0 3 ) 2 .4H 2 0, 0.8 mM MgS0 4 .7H 2 0, 0.4 mM 
KH 2 P0 4 , and 0.036 mM NaFe(iii).EDTA of macronutri- 
ents, and 18.50 uM H 3 B0 3 , 3.66 uM MnCl 2 .4H 2 0, 0.31 
uM ZnS0 4 .7H 2 O, 0.41 uM Na 2 Mo0 4 .2H 2 0, and 0.13 
uM CuS0 4 .5H 2 0 of micronutrients, at pH 6.0). One 32 



x 45 cm tub was sufficient to grow up to 200 plants. The 
nutrient solution was continually mixed with an aquar- 
ium pump and replaced weekly. The plants were grown 
in a glasshouse under natural light (12 hr day, 12 hr 
night) for 3 weeks and then treated with excess B by 
supplementing the aforementioned nutrient solution 
with 1 mM H 3 B0 3 . 

Quantification of boron tolerance 

Visual phenotyping was performed 3 weeks after B treat- 
ment when the three oldest trifoliate leaves (of about 
10-15 leaves) of each plant were scored using the following 
0-5 scale: 0 = green leaf, 1 = slight chlorosis, 2 > 10% leaf 
chlorosis, 3 = slight necrosis, 4 = 10-20% necrosis, 5 > 20% 
necrosis. The three leaves per plant were then harvested, 
pooled, and dried in a 65°C oven overnight. 

The dried leaf samples from parental lines and a sub- 
set of approximately 40 highly B toxicity tolerant and in- 
tolerant F 2 progeny from each cross were subjected to 
ICP-OES analysis [31] to measure their B concentra- 
tions. In an independent hydroponics experiment, root 
and leaf samples from the two tolerant and the intoler- 
ant parents, four plants each, were separately harvested 
in the manner as described above for ICP analysis. The 
wet roots were not rinsed but only blotted between 
paper towels to avoid B leaching as previously described 
[11]. The youngest two leaves were harvested to extract 
high quality genomic DNA for genotypic analysis. 

Quantification of plant biomass 

The weight of total plant dry matter [mg plant" 1 ] was 
measured from plants under high B stress in the txC F 2 
population using 137 F 2 lines and 5-9 repeats of their 
corresponding parents. Visual leaf symptom scores were 
taken before harvest, eight weeks after B treatment when 
the plants were a total of 11 weeks old. 

Development and detection of DNA markers 

DNA was extracted from freeze-dried leaf tissue using a 
previously documented protocol [32] with modifications 
as described [33]. To identify M. truncatula homologues 
of AtBORl and AtNIP5;l, their corresponding cDNA 
(GenBank: AB073713 and NM_1 17106) and protein 
(GenBank: BAC20173 and NP_192776) sequences were 
used to query the M. truncatula pseudomolecule data- 
base, version Mt3.5 (http://www.medicagohapmap.org). 
Once identified, nearby SSR markers were screened by 
PCR on a subset of lines from both txP and txC F 2 pop- 
ulations to establish linkage. PCRs with the SSR primers 
were carried out in 12.5 uL reactions containing IX re- 
action buffer (QIAGEN, Hilden, Germany), IX Q solu- 
tion (QIAGEN), 0.16 mM MgCl 2 , 0.32 uM of each 
primer, 0.4 U of Taq DNA polymerase (QIAGEN) and 
approximately 50 ng of DNA as template. The PCRs 
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were effected using the following cycling conditions: 90 
s at 95°C, then 35 cycles of 20 s at 95°C, 30 s at 58°C, 
and 30 s at 72°C, and a final extension step of 5 min at 
72°C All PCR products were electrophoresed on 8% 
polyacrylamide gels which were stained with ethidium 
bromide and visualised under UV light. 

For CAPS (cleaved amplified polymorphic sequence) 
marker development, PCR primers were designed on the 
sequence of the A17 genotype. Genomic regions were 
amplified by PCR in 12.5 \iL reactions as described earl- 
ier but with a 72°C extension step ranging from 30-90 s 
depending on the expected size of the amplicon. PCR 
products were resolved on 1.2% agarose gels and se- 
quenced using the Big Dye Terminator Version 3.1 label- 
ling kit (Applied Biosystems, Foster City, CA) according 
to the manufacturers protocol but using one-eighth of 
the recommended reaction mix. DNA sequencing was 
performed on an ABI 3730 DNA Sequencer (Applied 
Biosystems) and DNA sequences were analysed using 
Vector NTI software (Invitrogen). The SNP2CAPS pro- 
gram [34] was then used to develop CAPS markers from 
the sequence alignments. PCR products were digested with 
the relevant restriction enzymes according to the manufac- 
turer s protocol and separated on 1.5% agarose gels. 

Southern hybridization of parental DNA was carried 
out according to standard protocols [35] using 10 (ig of 
digested DNA per lane and the restriction enzymes 
BarnHI, Dral, EcoRI, EcoRV, Hindlll and Xbal (5 U/jig). 

Analysis of candidate B transporter gene expression 

To monitor the expression of candidate B transporter 
genes, plants were grown in a controlled- environment 
room for 3 weeks using the supported hydroponics system 
previously described [36] but with the same nutrient solu- 
tion as described earlier. Whole roots and the three youn- 
gest trifoliates from each plant were harvested and frozen 
in liquid nitrogen immediately before supplementing the 
nutrient solution with 1 mM H3BO3 (+B) and then at 3, 6, 
12, 24, and 48 hr after B treatment. A control set of plants 
was allowed to grow for 3 weeks after B treatment to en- 
sure that adequate symptoms developed. 

Total RNA was extracted from frozen root and young 
leaf tissue using TRI Reagent (Sigma- Aldrich, St. Louis, 
MO). For each time point, RNA was extracted from 
three independent biological replicates which were 
pooled to average out biological variability. To prepare 
cDNA template for quantitative real-time RT-PCR (Q- 
PCR), RNA was first treated with 2 U RNase-free DNase 
I (Ambion, Austin, TX) to remove residual DNA. cDNA 
was then synthesised from 3 \ig of RNA using 200 U of 
Superscript III reverse transcriptase (Invitrogen) 
according to the manufacturer s instructions with quality 
being confirmed via PCR detection of transcripts from 
the gene actin using earlier described primers [37]. Q- 



PCR was performed in an RG 3000 Rotor-Gene Real 
Time Thermal Cycler (Corbett Research, Sydney, 
Australia). The final reaction volume of 8 \iL comprised 
a mixture of IX Sensimix SYBR green PCR reagent 
(Bioline), 0.4 [iM of each forward and reverse gene- 
specific primer, and 3 \iL of a 1:20 dilution of cDNA. 
The PCR was effected using the following cycling condi- 
tions, 10 min at 95°C, followed by 35 cycles of 20 s at 
95°C, 30 s at 58°C and 30 s at 72°C. All Q-PCR reactions 
were carried out in duplicate and melt curve analysis 
was performed at the end of each run to confirm that 
there was no signal from non-specific binding products. 
No-template controls were also included in each run to 
test for possible contamination of assay reagents. The se- 
quences of all target gene primers and the reference 
gene actin primers are shown in additional file 3. 

To quantify gene expression levels, DNA standards were 
generated by cloning each target and reference gene PCR 
product into the pDrive cloning vector (QIAGEN) 
according to the manufacturers instructions. After spectro- 
photometric determination of plasmid DNA concentration, 
the vector was linearised with the restriction enzyme Nrul, 
and the copy number of standard DNA molecules calcu- 
lated using the formula (x g/\iL DNA / (plasmid length in 
base pairs x 660)) x 6.022 x 10 23 = y molecules/ (iL. Serial 
ten-fold dilutions of a given DNA standard (10 7 -10 1 copies/ 
\iL) were then incorporated into each run to quantify the 
number of target gene transcripts. These values were then 
normalised against the quantities of transcripts detected 
from the reference gene actin. 

Additional files 



Additional file 1: MtNIP3 cDNA sequence alignments for parental 
lines Caliph, Paraggio, and the intolerant tap. Single nucleotide 
polymorphisms between sequences are highlighted. 

Additional file 2: MtNIP3 protein sequence alignments for parental 
lines Caliph, Paraggio, and the intolerant tap. Amino acid differences 
between sequences are highlighted. 

Additional file 3: Q-PCR primers, primers for Southern probe and 
promoter analysis. 
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